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ABSTRACT 

We identify SDSS J010657.39- 100003.3 (hereafter J0106-1000) as the short- 
est period detached binary white dwarf (WD) system currently known. We targeted 
J0106— 1000 as part of our radial velocity program to search for companions around 
known extremely low-mass (ELM, ~ O.2M0) WDs using the 6.5m MMT. We detect 
peak-to-peak radial velocity variations of 740 km s^^ with an orbital period of 39.1 
min. The mass function and optical photometry rule out a main-sequence star compan- 
ion. Follow-up high-speed photometric observations obtained at the McDonald 2.1m 
telescope reveal ellipsoidal variations from the distorted primary but no eclipses. This 
is the first example of a tidally distorted WD. Modeling the lightcurve, we constrain 
the inclination angle of the system to be 67° ± 13°. J0106— 1000 contains a pair of 
WDs (0.17Mq primary -I- 0.43Mq invisible secondary) at a separation of O.32R0. The 
two WDs will merge in 37 Myr and most likely form a core He-burning single subd- 
warf star. J0106— 1000 is the shortest timescale merger system currently known. The 
gravitational wave strain from J0106— 1000 is at the detection limit of the Laser In- 
terferometer Space Antenna (LISA). However, accurate ephemeris and orbital period 
measurements may enable LISA to detect J0106— 1000 above the Galactic background 
noise. 

Key words: (stars:) binaries (including multiple): close — (stars:) white dwarfs — 
(stars:) individual (SDSS J010657.39- 100003.3) — Galaxy: stellar content 



1 INTRODUCTION 

ELM WDs are ideal targets for finding binary WD merger 
systems. Short period binary stars interact early in their 
stellar evolution, experience enhanced mass-loss during one 
or two common-envelope phases ( Sarna et al. .l996i ) , and end 
up as lower mass WDs. Thus the most compact binary sys- 
tems are expected to form ELM WDs, and a survey targeting 

EL M WDs should discover mergin g sys tems. 

iKilic et all (|2009l . I2OI0I . I2OIII ) and iBrown et al.) (|2010l ') 
have established a radial velocity program, the ELM Sur- 
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vey, to search for companions around known ELM WDs in 
the SDSS Data Release 7 footprint. The discovery of 12 bi- 
nary WD merger systems in a sample of two dozen WDs 
observed to date has tripled the number of known binary 
WD merger systems. Eighteen of the M ^ 0.25Mq WDs in 
the ELM Survey are in 1-24 hr period binaries with merger 
times as short as 100 Myr. Six of these systems have ex- 
treme mass ratios {M1/M2 = q ~ 0.2), which may lead to 
stable mass transfer AM CVn systems. If the mass-accreting 
WDs in these systems are massive, they can potentially form 
Type la supernovae (SNe, IWebbinklflQSi : llben fc Tutukovl 
|T984). Alternatively, accretion of helium from a companion 
may lead to the detonation of the surface helium layer on 
a C/0 WD in a fas t and faint supernova, i.e. SNe ".la" 
(|Bildsten et allbOOTl L 

Here we present the exciting discovery of a new binary 
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system found in the ELM survey. J0106— 1000 was originally 
classified as a sub dwarf s tar in the SPSS D R4 WD catalog of 
lEisenstein et alj ()2006l ) . iKleinmanI ()2010l ) re-classified it as 
an ELM WD based on a reanalysis of the SDSS spectroscopy. 
Our follow-up radial velocity and high speed photometric 
observations demonstrate that J0106— 1000 contains a pair 
of WDs with an orbital period of only 39.1 minutes. This sys- 
tem, the shortest period detached binary WD system known, 
presents the first detection of a tidally distorted WD. 

In Section 2 we describe our spectroscopic and photo- 
metric observations. In Sections 3 and 4 we constrain the 
physical parameters of the binary and discuss the nature 
and future evolution of the J0106— 1000 system. We con- 
clude in Section 5. 



2 OBSERVATIONS 

We used the 6.5m MMT with the Blue Channel spec- 
trograph to obtain medium resolution spectroscopy of 
J0106-1000 on UT 2010 Dec 1 - 3. We operate the spec- 
trograph with the 832 line mm~^ grating in second order, 
providing wavelength coverage from 3600 A to 4500 A and 
a spectral resolution of 1.2 A. We obtain all observations 
at the parallactic angle, with a comparison lamp exposure 
paired with every observation. We flux-calibr ate using blue 
spectrophotometric standards (Mas sev et al.| [l988). and we 
measure radial velocities using the cross-correlation package 
RVSAO. Th e details of our data re duction procedures are 
discussed in iKilic et all (|2009l . I2OI0I ). We check the stabil- 
ity of the spectrograph using the Hg line at 4358.34 A from 
Tucson/Nogales street lights. We measure an average veloc- 
ity offset of —0.9 ± 0.3 km s^^ for this line over the entire 
three nights of observations. 

J0106-1000 is relatively faint {g = 19.8 mag). We 
started our observations with 10 min exposures. After de- 
tecting > 400 km s~^ velocity variations in 15 min, we de- 
creased the individual exposure times to 8 min. Realizing 
that J0106— 1000 is a very short period system after our 
second night of observing, we also acquired high speed pho- 
tometric observations of J0106— 1000 using the McDonald 
2.1m Ot to Struve Telescope with the Argos frame transfer 
camera (|Mukadam fc Natheill2005l ) on UT 2010 Dec 3. Ar- 
gos provides a field of view of 2.8' x 2.8'. We obtained time 
series photometry of J0106— 1000 with the BG40 filter every 
30 s for about 2.6 hr. The Argos field-of-view includes several 
comparison stars that are useful for relative photometry. 



3 RESULTS 

3.1 The Orbital Period 

Table 1 lists our radial velocity measurements for 
J0106 — 1000. We compute best-fi t orbital elements using the 
code of lKenvon fc Garcia 1 1| 19861 ). which weights each veloc- 
ity measurement by its associated error. The uncertainties in 
the orbital elements are derived from the covariance matrix 
and X^- To verify these uncertainty estimates, we perform 
a Monte Carlo analysis using 10'' sets of simulated radial 
velocities. We adopt the inter-quartile range in the period 
and orbital elements as the uncertainty. 

J0106— 1000 exhibits radial velocity variations with a 



Table 1. Radial Velocity Measurements for J0106-1000 
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Figure 1. The radial velocities of J0106— 1000 observed over 
three nights in 2010 December (left panels). The right panel shows 
all of these data points phased with the best-fit period. The solid 
line represents the best-fit model for a circular orbit with a period 
of 39.1 min and K = 369.5 km s~^. 

semi- amplitude of K = 369.5 ± 3.6 km s~'and orbital pe- 
riod of P = 0.027153 ± 0.0000195 d, or 39.100 ± 0.028 min. 
Figure 1 shows the best-fit orbit compared to the observed 
radial velocities. The relatively long exposure times (8 min) 
compared to the orbital period (39.1 min) results in an un- 
derestimated velocity semi-amplitude, K. This is a direct 
consequence of the sine curve not being linear when the ve- 
locities are at the extremes. To verify this effect, we sampled 
a sine curve at the exact 21 phases of our observations with 
P/5 long integrations. We recover the exact period, but K is 
systematically underestimated by 6.5%. Thus, the corrected 
velocity semi-amplitude for J0106— 1000 is K — 395.2 km 
s~^. With this correction, J0106— 1000 has a mass function 
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Figure 2. Model fits (red lines) to the Balmer line profiles of 
J0106— 1000 (jagged lines, top panel). The spectral energy distri- 
bution of J0106-1000 (based on the SDSS photometry, filled cir- 
cles) compared to the best- fit model (open circles, bottom panel). 



of / = 0.1736 ± 0.0047 M0. The systemic velocity (after 
subtracting the gravitational redshift of 1.9 km s~^) is 0.3 
± 2.7 km s~^ and the time of spectroscopic conjunction is 
HJD 2455531.633574 ± 0.000129 d. 



3.2 The Physical Parameters of the Binary 

Our time-series spectroscopy provides for robust determina- 
tions of effective temperature and surface gravity. We per- 
form stellar atmosphere model fits using synthetic DA WD 
spectra kindly provided by D. Koester. The grid of WD 
model atmospheres covers effective temperatures from 6000 
K to 30,000 K in steps of 500 K to 2000 K, and surface grav- 
ities from logg = 5.0 to 9.0 in steps of 0.25 dex. We perform 
fits to the Balmer line profiles using the average composite 
spectra. We also perform fits to the individual spectra to de- 
rive a robust statistical error estimate. Figure 2 shows the 
observed Balmer line profiles (jagged lines) for J0106— 1000 
compared to our best-fit model (solid line). The best-fit 
model has T^a = 16485 ± 456 K and logg = 6.01 ± 0.04. 

Previous stellar atm osphere fits find a comparable log g 
but a 2500 K hotter T^s l|Kleinmanll2010l ). The discrepancy 
in their temperature estimate is most likely due to the low 
signal-to-noise SDSS spectrum. Our best-fit model matches 
our S/N = 40 spectrum and the SDSS photometry (Figure 

2). r— — 1 

Based on th e improved IPanei et al.l ()2007[ ) tracks (see 
iKilic et al.|[201Gl ) for ELM WDs, J0106-1000 is a 1.1 Gyr 
ol(£] O.17M0 WD with a radius R = 0.057 Rq. Its ab- 
solute magnitude Mg = 7.8 corresponds to a distance 



^ This age estimate is somewhat uncertain due to the assumption 
on the thickness of the surface hydrogen layer for 0.17MqWDs. 



of 2.4 kpc . Based on fi ve ep ochs from the USNO-B and 
the SDSS, iMunn et al.l (|200^ ) measure a proper motion of 
(paCosS,fj.s) = (20.2,-10.5) mas yr"^ J0106-1000 is 2.3 
kpc below the Galactic plane and it has U = — 115±43,V — 
-222±43, and W = -15±12 km s~^ with respect to the lo- 
cal standard of rest (|Hogg et al.ll2005l ). Clearly, J0106-1000 
is a halo star. 

The mass function implies a companion mass of ^0.37 
Mq. The orbit is far too small to contain a main-sequence 
star of 0.37 Mq or more . Therefore, the companion is 
a compact object. Based on the mass function alone, the 
probability of a neutron star (1.4-3.0 Mq) companion is 
7.3%. The probability of a SNe la, for which the companion 
would be a 1.23 - 1.40 Mq WD, is only 1.9%. However, if 
sub-Chandrasekhar mass WDs do explode as Type la SNe 
(Ivan Kerkwiik et"aLll2010l ). this probability may be higher. 



3.3 The Light Curve 

The chance of an eclipse is relatively high for J0106— 1000. 
For an edge-on orbit, the companion would be a 0.37Mq 
WD. To avoid detection in the SDSS photometry, we assume 
that the companion is 10 x fainter than the visible WD and 
thus it has T^s 14600 K and _R 0.02Rq. The eclipse 
depth and duration would be 12% and ~ 90 s (assuming a 
total eclipse), respectively. Due to the relatively large size 
of the visible WD compared to the orbital separation, the 
probability of a grazing eclipse is 25%. 

Figure 3 shows the Argos light curve of J0106— 1000 
(top panel) over four orbits. J0106— 1000 is relatively faint 
and the Argos light curve has a few percent scatter. We 
do not detect any pulsations at the ^ 0.8% level. However, 
the fourier transform of the J0106— 1000 light curve reveals a 
significant peak at 1187 ± 13 s (half the orbital period) with 
an amplitude of 1.7% ± 0.3%. We use the ephemeris from 
the radial velocity observations obtained on the same night 
to calculate the phase for our photometric observations. The 
middle panel in Figure 3 shows the light curve folded over 
the best-fit orbital period. There are essentially four separate 
observations every 30 s in the folded light curve. The bottom 
panel in Figure 3 shows the same light curve binned by four 
points. This panel clearly shows the ellipsoidal variations; 
the ELM WD is distorted due to the companion. This is the 
first detection of ellipsoidal variations for a WD. 

To verify that the observed variations in the light curve 
are not a statistical fiuctuation, we perform a bootstrap 
analysis. We randomly permute the light curve points in 
time and create 10^ simulated light curves. This analysis 
shows that the probability of the measured ellipsoidal vari- 
ation signal being a random one is smaller than 10~^. We 
also checked for variations in the brightest reference star by 
folding its light curve for the 39.1 min period of J0106— 1000. 
This analysis did not reveal any periodic variations in the 
reference star. These tests show that the detected variability 
is highly significant. 

To model th e ellipsoidal variati o ns an d the reflection 
effect, we follow i Morris fc NaftilanI l| 19931 ). In their for- 
malism, the amplitude of the ellipsoidal effect is roughly 
Sfeii = (m2/mi)(ri/a)^, where a is the orbital semi-major 
axis and ri is the ra dius of the primary (jZucker et al. I l2007l : 
IShporer et a"l]|2010l ). For J0106-1000, (m2/mi) « 2.2 for 
an edge-on orbit, yielding S/m = 1.4%. For a more accurate 
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Figure 3. High speed photometry of J0106— 1000 over four or- 
bital periods (top panel). The middle and bottom panels show 
the same light curve folded over the orbital period and binned 
by four points, respectively. Our best- fit model (i = 73°) and a 
comparable model (i = 60°) are shown as solid and dotted lines. 



estimate, we follow equation 1 in iMorris fc NaftilanI (|l993l ) 
and include the low order terms up to the cos 4(/) term, 
where (j) = 0° when the primary star is farthest from the 
observer. The transmission curve of the Argos camera with 
the BG40 filter is similar to a B-band filter. We use limb 
darkening and gravity-darkening coefficients of ui = 0.36 
and Ti = 0.487, respect ively (see equation 1 and Table 1 in 
iMorris fc Naftilanll 19931 ). The results change only slightly for 
different limb-darkening and gravity-darkening coefficients. 
The ellipsoidal variations are dominated by the cos 2(j} term, 
which has an amplitude of 1.8% for the J0106— 1000 system 
at 90° inclination. We estimate the contribution from the 
reflection effect to be <0.1% and therefore neglect it. On 
the other hand, t he relativistic beam i ng effect is expecte d 
to be about 0.3% (|Maxted et al.ll2000l : [Shporer et al.ll20ld ). 
We include the beaming effect, but not the eclipses in our 
calculations. 

The amplitudes of the ellipsoidal variations and the rel- 
ativistic beaming effect depend on the companion mass and 
orbital separation, which depend on the inclination of the 
system. Hence, the predicted amplitudes can be represented 
by a single parameter, the inclination angle. We create model 
light curves for each inclination angle and chose the model 
with the least a-s the best-fit model (shown as a solid line 
in Figure 3). This model has an inclination angle of i = 73° 
and it matches the phase and amplitude of the variations rel- 
atively well. The absence of eclipses in the photometry is also 
consistent with this inclination angle. The dotted line shows 
a comparable model with an inclination angle of i = 60°. 
The similarities between the two models and the relatively 
large scatter in our photometry indicate that the error in 
inclination is relatively large. 



To constrain the inclination angle of the system more 
accurately, we perform a Monte Carlo analysis where we 
replace the measured flux / with f + g 5f, where 5f is the 
error in flux and p is a Gaussian deviate with zero mean 
and unit variance. For each of the 10^ sets of simulated light 
curves, we repeat our analysis and find the best-fitting model 
(and the inclination angle). We adopt the interquartile range 
as the uncertainty. The Monte Carlo analysis shows that 
the J0106— 1000 light curve is best explained by a model 
with i — 67° ± 13° . Therefore, the companion is most likely 
a 0.43Mq object at an orbital separation of 0.32Rq. The 
mass ratio of the binary is q = 0.4 and the merger time due 
to gravitational wave radiation is 37 Myr. 



4 DISCUSSION 

J0106— 1000 is the shortest period detached binary WD 
system currently known. Even though its orbital period is 
comparable to the AM CVn systems, its spectrum shows 
only hydrogen absorption lines. There is no evidence of in- 
teraction or mass accretion between the two components 
other than the slight distortion of the ELM WD due to mu- 
tual gravitation. The observed ellipsoidal variations are ex- 
tremely useful for constraining the inclination angle and the 
masses of both components of the system. The optical spec- 
troscopy and photometry of this system is best explained 
by a binary system containing two He-core WDs with Mi — 
0.17Mq and Ah = 0.43Mq at a separation of O.32R0. 

The future evolution of the system depends on the 
mass ratio of the two components. For a mass ratio of 
q = 0.4, J0 106— 1000 will lik ely have unstable mass transfer 
and merge. iDan et al.l l|201lf ) simulate the mergers of double 
degenerate systems including O.2M0 ELM WDs with 0.3- 
0.8Mq WD companions. Their smoothed-particle hydrody- 
namic models indicate that systems with q ^ 0.25 have un- 
sta ble mass transfer and merge after 20 -80 orbits (although 
see lMotl et al.ll2007l : iRacine et al.ll2007l , for the possibility of 
stable mass transfer). Unstable mass transfer may also lead 
to the detonation of the surface hel ium layer on a C/ Q WD 
via Kelvin-Helmholtz instabilities. ICuillochon et af] (|2010l ) 
find that the required conditions for triggering a surface ex- 
plosion are only achieved in ^ 0.8Mq WD accretors with 
>O.2M0 companions. Hence, no efficient carbon burning is 
expected in binary WD systems containing 0.45-0.8Mq C/O 
WD accretors and ELM WD donors. Based on these results, 
J0106— 1000 will likely merge and create a 0.6Mq core- He 
burning subdwarf in 37 Myr. This mass is close t o the canon- 
ical mass of O.5M0 for subdwarfs (|Hebeij|2009l ). 

iBrown et al.l (|201ll ) estimate the merger rate of ELM 
WDs from a complete, color-selected sa mple of ELM WDs 
found in the Hypervelocity Star Survey (iBrown et al.|[2006l . 
l2009t ). Roughly 5 xlO** (with a factor of few uncertainty) 
ELM WDs formed in the Galactic disk in the last Gyr. 
About 70% of these systems merge in less than a Gyr, and 
66% of these systems have mass ratios q ^0.25. From these 
estimates, roughly 2.3 x 10"* ELM WDs merged in the past 
Gyr. There are three stars with q ^0.25 and merger times 
shorter than a Hubble time in their sample. These three 
stars, J0818-I-3536, J0923-I-3028, and J1053-I-5200, contain 
O.17-O.23M0 -I- O.33-O.44M0 WDs (assuming an average 
inclination angle of 60°). Based on the evolutionary calcu- 
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lations by iDan et al] l|201 j ). they are likely to form 0.5- 
0.67Mq single subdwarfs. Hence, the formation rate of sin- 
gle subdwarfs through mergers of ELM WDs is roughly 2.3 
xlO* in the last Gyr. Adding J0106-1000 to this sample 
would increase this rate by about 50% due to its relatively 
sho rt merger t i me. 

iNelemand (|2010l ) presents population synthesis models 
for the Galactic population of subdwarf B stars. He pre- 
dicts a total number of 5.6 x 10^ single subdwarfs to form 
as a result of He WD mergers. The resulting mass distri- 
bution is centered around 0.5Mq with a tail toward higher 
masses, similar to the mass distribution of the merger sys- 
tems discussed above. Hence, the ELM WD merger systems 
contribute significantly to the population of single subdwarfs 
in the Galaxy. 

Short peri od binary WDs a re imp ortant gravitational 
wave sources. iNelemans et al.l ((200J) and [Roclofs ct al.l 
l|2007l l argue that about half a dozen AM CVn binaries 
should be detected by LISA. With an orbital period sim- 
ilar to the known AM CVn systems, J0106— 1000 may 
be a promising candidate for detection. The orbital pe- 
riod, inclination, and model-dependent distance estimate 
for J0106— 1000 yield the gravitational wave strain at 
Earth, log ft = —22 .7 at a frequency logz/ (Hz) — —3.07 
iRoelofs et al.ll2007l 'l. This is at the S/N = 1 detection limit 
of LISA after 1 year of observations. Confusion with Galactic 
noise sources decreases with accurately known orbital peri- 
ods. Our ephemeris and orbital period measurements may 
enable LISA to detect J0106— 1000 above the Galactic back- 
ground noise after a few years of observations. 



5 CONCLUSIONS 

We discovered the shortest period detached binary WD sys- 
tem currently known. This system also presents the first 
detection of a tidally distorted WD. We constrain the in- 
clination angle of the system using high-speed photometric 
observations. J0106— 1000 contains a pair of low-mass WDs 
at an inclination angle of 67° ± 13°. Follow-up high-speed 
photometric observations at a larger telescope will be useful 
to better constrain the inclination (and therefore the com- 
panion mass) and to search for grazing eclipses. The two 
WDs will merge in 37 Myr and most likely form a core He- 
burning single subdwarf star. 
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